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ON  THE  INITIAL  STAGES  OF  VORTEX  WAVE  INTERACTIONS 
IN  HIGHLY  CURVED  BOUNDARY  LAYER  FLOWS 


Philip  Hal P 

Department  of  Mathematics 
University  of  Manchester 
Manchester  M13  9PL 
UNITED  KINGDOM 

ABSTRACT 

The  nonlinear  interaction  equations  describing  vortex- Rayleigh  wave  interactions  in  highly 
curved  boundary  layers  are  derived.  These  equations  describe  a  strongly  nonlinear  interac¬ 
tion  between  an  inviscid  wave  system  and  a  streamwise  vortex.  The  coupling  between  the 
two  structures  is  quite  different  than  that  found  by  Hall  and  Smith  (1991)  in  the  absence 
of  wall  curvature.  Here  the  vortex  is  forced  over  a  finite  region  of  the  flow  rather  than  in 
the  critical  layer  associated  with  the  wave  system.  When  the  interaction  takes  place  the 
wave  system  remains  locally  neutral  as  it  moves  downstream  and  it’s  self  interaction  drives 
a  vortex  field  of  the  same  magnitude  as  that  driven  by  th~  ,vall  curvature.  This  modification 
of  the  mean  state  then  alters  the  wave  properties  and  forces  the  wave  amplitude  to  adjust 
itself  in  order  that  the  wave  frequency  is  constant.  Solutions  of  the  interaction  equations 
are  found  for  the  initial  stages  of  the  interaction  in  the  case  when  the  wave  amplitude  is 
initially  small.  Our  analysis  suggests  that  finite  amplitude  disturbances  can  only  exist  when 
the  vortex  field  is  finite  at  the  initial  position  where  the  interaction  is  stimulated. 
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1.  Introduction 


Our  concern  is  with  the  interaction  of  Rayleigh  waves  and  streamwise  vortices  in  pressure 
gradient  driven  boundary  layers  on  highly  curved  walls.  The  first  step  is  to  derive  the  coupled 
strongly  nonlinear  interaction  equations  governing  the  small  wavelength  vortex  field  driven 
by  small  amplitude  Rayleigh  waves.  The  interaction  equations  we  obtain  differ  from  those 
found  previously  by  Hall  and  Smith  (1991).  In  the  latter  paper  it  was  found  that  three- 
dimensional  Rayleigh  waves  drive  the  vortex  field  in  the  critical  layer.  The  vortex  field  then 
acts  back  on  the  wave  indirectly  because  the  wave  satisfies  a  modified  form  of  Rayleighs 
equation  which  depends  on  the  spamvise  variation  of  the  mean  flow.  Here  we  find  that 
the  Rayleigh  waves  cannot  drive  the  vortex  field  in  the  critical  layer.  Instead  we  find  that 
the  forcing  is  distributed  over  the  region  where  the  vortex  activity  occurs,  this  enables  us 
to  reduce  the  interaction  problem  to  a  nonlinear  partial  differential  system  independent  of 
the  spanwise  variable.  We  construct  the  small  amplitude  solutions  of  this  system  in  two 
situations  and  show  that  the  possible  emergence  of  the  vortex-wave  interaction  is  sensitively 
controlled  by  the  nature  of  the  incoming  mean  profile. 

The  different  stages  of  vortex-wave  interactions  associated  with  viscous  (Tollmien- 
Schlichting)  travelling  waves  are  described  in  the  papers  by  Hall  and  Smith  (1988,  1989, 
1990,  1991),  Smith  and  Walton  (1989),  Bassom  and  Hall  (1989),  Bennett,  Hall  and  Smith 
(1991)  and  Blackaby  (1991).  Typically  it  is  found  that  Tollmien-Schlichting  waves  force  the 
vortex  field  in  a  thin  viscous  layer  at  the  wall  for  external  flows  and  away  from  the  walls  for 
fully  developed  internal  flows.  For  external  flows  the  forcing  results  in  an  inhomogeneous 
boundary  condition  for  the  vortex  velocity  component  in  the  spanwise  direction.  The  sit¬ 
uation  with  Rayleigh  waves  is  more  complicated  and  Hall  and  Smith  (1991)  showed  that 
Rayleigh  waves  force  the  vortex  in  the  critical  layer  associated  with  the  wave.  More  pre¬ 
cisely  the  wave  induces  a  finite  jump  in  the  gradient  of  the  spanwise  velocity  component  of 
the  vortex  across  the  critical  layer.  The  solution  of  the  interaction  equations  in  this  case 
is  made  difficult  by  the  spanwise  and  streamwise  variation  of  the  critical  layer  position;  as 
yet  no  numerical  solutions  of  the  interaction  equations  derived  by  Hall  and  Smith  have  been 
obtained. 

In  the  present  paper  we  examine  the  Rayleigh  wave-vortex  interaction  problem  in  the 
case  of  highly  curved  boundary  layers.  This  regime  corresponds  to  high  values  of  the  Gortler 
number  associated  with  the  flow  and  we  assume  that  the  vortex  wavelength  is  small  compared 
to  the  boundary  layer  thickness.  The  latter  assumption  enables  us  to  make  use  of  the  small 
wavelength  large  amplitude  structure  of  vortices  discussed  by  Hall  and  Lakin  (1988).  The 
latter  calculation  was  a  development  of  the  linear  and  weakly  nonlinear  investigations  of 
the  Gortler  vortex  mechanism  by  Hall  ( 1982a, b).  In  the  present  calculation  we  show  that 
at  small  wavelengths  the  critical  layer  forcing  of  the  vortex  by  Rayleigh  waves  becomes 
exponentially  small  compared  to  a  new  distributed  mechanism  operational  over  a  finite  part 
of  the  boundary  layer.  We  shall  see  that  the  forcing  is  confined  to  the  finite  part  of  the 
boundary  layer  where  vortex  activity  occurs. 

Th_'  approach  of  Hall  and  Lakin  (1988)  applied  to  the  present  problem  shows  that  the 
interaction  equations  can  be  significantly  reduced  to  a  stage  where  the  spanwise  dependence 
of  the  disturbances  is  described  analytically.  The  means  that  the  strongly  nonlinear  interac¬ 
tion  equations  for  vortex-wave  interactions  at  small  wavelengths  can  be  reduced  to  a  partial 
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differential  system  dependent  only  on  the  two-dimensional  boundary  layer  variables.  These 
interaction  equations  are  derived  in  §2. 

A  non-trivial  problem  associated  with  the  equations  governing  vortex-wave  interactions 
is  the  determination  of  appropriate  initial  conditions  for  the  equations.  In  §3  we  tackle 
the  problem  for  the  case  when  a  Rayleigh  wave  of  small  amplitude  is  generated  from  a 
finite  amplitude  vortex  field.  We  shall  derive  constraints  on  the  vortex  field  which  allow 
for  the  emergence  of  a  vortex- wave  interaction.  In  §4  we  investigate  the  possibility  of  the 
spontaneous  generation  of  Rayleigh  waves  and  vortices  at  some  position  in  the  streamwise 
direction.  Our  results  suggest  that  the  stringent  constraints  which  must  be  satisfied  if  such 
an  interaction  is  to  take  place  virtually  rule  out  the  spontaneous  generation  of  Rayleigh 
waves  and  vortices.  In  §5  we  discuss  the  results  of  §3,4  and  draw  some  conclusions. 

2.  Derivation  of  the  Vortex- Wave  Interaction  Equations  for  Highly  Curved  Flows 

We  consider  the  flow  of  a  viscous  incompressible  fluid  of  density  p,  kinematic  viscos¬ 
ity  v  past  a  rigid  wall  defined  by  y  =  0  with  respect  to  a  coordinate  system  ( x ,  y,  z )  with 
x,y,z  measuring  distance  along  the  wall,  normal  to  the  wall  and  in  the  spanwise  direc¬ 
tion  respectively.  We  assume  that  x,y,z  have  been  made  dimensionless  using  the  lengths 
L,  LR -1/2,  Li?-1/2  where  L  is  a  lengthscale  in  the  x-direction  and  R  is  defined  by 

R=—.  (2.1) 

v 

Here  U0  is  a  typical  fluid  speed  in  the  x  direction  and  the  Reynolds  number  R  is  taken  to 
be  large.  If  we  take  the  wall  curvature  to  be  a-1x(x)  then  the  equations  to  determine  the 
combined  Gortler  vortex-mean  flow  are 

V  •  w*  =  0, 

l  0  \  (  -px\  (2.2) 

w*  •  Vm*  +  xGu*2/ 2  =  -pi  +  Aw*. 

V  0  1  \  -PW 

In  the  above  equations  G  is  the  Gortler  number  defined  by 

G  =  R''2j, 

f  j2 - 

which  is  held  fixed  in  the  limit  R  — ►  oo,  p--^—  is  the  streamwise  pressure  gradient  at  the 
edge  of  the  boundary  layer  and  u*,u*,u>*,p*  have  been  made  dimensionless  using  the  scales 
(/o,  (/oR-1^2,  UqR~1^2  and  p(J$R~^  respectively.  Finally  the  operator  A  appearing  in  (2.2) 
is  the  two-dimensional  Laplacian  d2  +  d].  Thus  streamwise  diffusion  is  negligible  in  (2.2) 
since  it  operates  on  a  longer  lengthscale  than  diffusion  in  the  t/,2  directions.  In  addition, 
p *  is  negligible  so  that  (2.2)  is  parabolic  in  x ,  therefore  no  mechanism  to  produce  upstream 
influence  is  present. 

In  order  to  study  the  influence  of  vortex  flows  on  the  inviscid  instability  of  shear  flows 
Hall  and  Horseman  (1991)  superimposed  on  the  flow  an  infinitesimally  small  Rayleigh  wave. 


Following  these  authors  we  therefore  replace  the  functions  u*,v*,w *  used  in  the  derivation 
of  (2.2)  by  writing 

u*  =  u(x,y,z)  +  £j75-  [U(x,y,z)E  +  complex  conjugate]  H - , 

v *  =  u(x,  y,  z)  -f  8  [V(x,  y,  z)E  -f  complex  conjugate]  +  ,  (2.3) 

w*  =  iu(x,  y,  z)  +  8  [W(x,  y,  z)E  -f  complex  conjugate]  +  ••■ , 

where 

E  =  exp  iR1 12  |  J  a(x)dx  —  fit  j  .  (2.4) 

Thus  we  have  assumed  that  the  disturbance,  of  arbitrarily  small  size  6 ,  is  periodic  in  time 
(with  t  scaled  on  R~1^2UqX)  and  varies  on  a  short  0(R~ 1//2)  lengthscale  in  the  x  direction. 
Note  that  if  8  is  sufficiently  small  then  u,v,w  satisfy  (2.2)  with  the  asterisks  removed.  The 
equations  satisfied  by  the  disturbance  in  the  limit  R  — *  oo  are  found  to  be 

iaU  +  Vy  +  Wt  =  0, 


ia(u  —  c)U  +  Vuv  +  Wuz  —  —iaP, 
ia(u  —  c)V  =  — Pv, 


(2.5) 


ia(u  —  c)W  =  —Pz, 


with  c  =  c(x)  =  fl/a  and  P  denoting  the  pressure  perturbation  corresponding  to  U,V,W. 
A  more  convenient  form  of  (2.5)  is  obtained  by  eliminating  the  velocity  field  U,  V,  W  to  give 
the  pressure  form  of  the  Rayleigh  equation  for  longitudinal  vortex  flows;  we  thus  obtain 

((^^)>  +  ((^ji),“(^j5  =  0,  (2'6a) 

and  this  must  be  solved  with  P  periodic  in  2  and  subject  to 


Py  =  0,  y  =  0, 

P  — *  0,  y  — ►  00. 


(2.66,  c) 


For  a  given  velocity  field  it(x,  y,  2)  the  eigenvalue  problem  (2.6)  can  be  solved  for  a  =  a(ft,  c), 
and  the  flow  is  unstable  if  eigenvalues  can  be  found  with  a*  <  0.  A  similar  eigenvalue  problem 
for  the  temporal  instability  problem  with  a  real  and  c  complex  can  also  be  used  to  classify 
a  given  flow  as  stable  or  unstable.  Using  the  latter  approach  Hall  and  Horseman  (1991) 
were  able  to  show  that  Gortler  vortices  cause  a  Blasius  boundary  layer  to  become  inviscidly 
unstable  when  they  become  of  sufficiently  large  amplitude.  The  results  found  by  the  latter 
authors  were  in  excellent  agreement  with  the  experimental  observations  of  Swearingen  and 
Black  welder  (1987). 

In  the  context  of  vortex-wave  interaction  theory  interest  centres  on  the  neutral  solutions 
of  (2.6)  which  have  a  =  a(x),  ac  =  fi  =  constant  and  8  =  6(x)  chosen  to  be  large  enough  to 
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maintain  the  wave  in  this  neutral  state  as  it  moves  downstream.  In  Hall  and  Smith  (1991), 
hereafter  referred  to  as  HS,  it  was  shown  that  the  forcing  takes  place  in  the  critical  layer  and 
leads  to  a  jump  in  wy  across  the  layer.  However  the  expression  for  the  jump  in  wy  vanishes 
if  P  is  independent  of  2  at  the  critical  layer.  This  would  for  example  be  the  case  in  the 
degenerate  case  u  =  u(x,y)  which  applies  when  no  vortices  are  present  in  the  flow.  Thus  a 
fundamental  property  of  the  interaction  is  that  it  can  only  occur  when  the  streamwise  flow 
has  a  spanwise  dependence. 

In  fact  there  exists  a  class  of  flows  which  support  large  amplitude  streamwise  vortices 
over  only  a  finite  part  of  a  boundary  layer.  If  the  critical  layer  associated  with  these  flows 
is  in  a  region  where  no  vortex  activity  occurs  then  the  mechanism  described  in  HS  cannot 
operate.  The  class  of  flows  referred  to  above  corresponds  to  the  strongly  nonlinear  Gortler 
vortex  flows  first  discussed  by  Hall  and  Lakin  (1988).  These  flows  correspond  to  the  limit 
1  (or  equivalently  a:>l)  and  have  the  structure  sketched  in  Figure  (1.1). 

In  Regions  I  and  II  there  is  no  vortex  activity  and  at  zeroth  order  the  mean  state  simply 
satisfies  the  equations  governing  the  unperturbed  boundary  layer.  In  Region  III,  which 
corresponds  to  yj  <  y  <  y2,  a  finite  amplitude  vortex  exists  and  drives  the  mean  state.  The 
vortex  activity  is  reduced  exponentially  to  zero  in  thin  shear  layers  surrounding  j/i  and  y2. 
It  follows  that  if  we  have  a  situation  in  which  the  critical  layer  of  the  system  (2.6)  associated 
with  the  above  flow  is  below  y  =  J/i,  then  an  alternative  mechanism  to  that  proposed  by 
HS  must  be  found.  For  the  moment  we  assume  that  u(x,y,z)  appearing  in  (2.2),  with 
(u*,v*,wm  ,p*)  replaced  by  (u,v,u>,  p)  has  a  region  adjacent  to  y  =  0  where  it  is  independent 
of  2  and  that  the  critical  layer  occurs  in  this  region.  We  now  have  a  situation  where  the 
Rayleigh  wave  satisfying  (2.6)  cannot  force  the  vortex  flow  in  the  manner  discussed  in  HS 
in  the  region  of  vortex  activity. 

In  fact  the  vortex-wave  interaction  in  such  a  case  occurs  in  a  simpler  manner  than  that 
found  in  HS.  In  order  to  allow  the  Rayleigh  wave  to  drive  the  vortex  we  let  8  in  (2.3) 
increase  until  the  nonlinear  terms  involving  {/,  V,  YV  in  the  y  and  2  momentum  equations 
become  comparable  with  those  present  for  the  vortex  field.  It  is  easy  to  see  that  this  occurs 
when  8  =  R~'!28  with  8  of  order  R°.  In  this  case  the  momentum  equations  in  (2.2)  are 
modified  to  give 


u  ■  Vu  -f 


0 


xGu2/ 2 
6 


l  0  \ 


~Vr  \ 
~Py 

+  -  82 

-ioUV  +  VVy  +  WVZ  +  c.c. 

-Pz  ) 

v  -iallW  +  VWy  +  WW:  +  c.c.  ) 

where  ‘c.c.1  denotes  ‘complex  conjugate.’ 

Thus  the  vortex  field  is  now  driven  by  the  Rayleigh  wave  through  the  forcing  terms  on 
the  right  hand  side  of  (2.7).  The  Rayleigh  wave  is  determined  by  (2.6)  and  so  is  itself  a 
function  of  u(x,y,z),  thus  there  is  a  strongly  nonlinear  coupling  between  the  vortex  and 
the  wave  fields.  In  order  that  the  Rayleigh  wave  remains  neutral  as  it  moves  downstream 
the  wavespeed  adjust  itself  to  the  fluid  speed  at  the  inflection  point.  If  the  frequency  is  to 
remain  fixed  then  there  is  an  apparent  difficulty  because  the  wavespeed  and  wavenumber 
are  effectively  fixed  by  the  previous  consideration.  The  required  extra  degree  of  freedom  is 
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found  by  allowing  S  to  adjust  itself  so  as  to  modify  the  mean  flow  in  such  a  way  that  the 
frequency  remains  fixed. 

We  note  that  (2.7)  fails  if  u  =  c  in  the  region  of  vortex  activity  and  in  this  case  we  return 
to  the  situation  discussed  by  HS.  A  further  important  point  to  notice  here  is  that  the  size  of 
the  Rayleigh  wave,  0{R~ 1/,2))  required  to  drive  the  interaction  is  so  small  that  the  Rayleigh 
wave  critical  layer,  of  depth  0(R~X /6),  remains  linear. 

We  now  indicate  how  the  short  wavelength  asymptotic  structure  given  by  Hall  and  Lakin 
(1988)  can  be  modified  to  take  account  of  the  forcing  terms  arising  from  the  interaction  of 
the  Rayleigh  wave  with  itself.  We  assume  that  the  vortex  activity  is  confined  to  the  region 
y\  <  y  <  y 2.  Here  we  take  the  vortex  wavenumber  to  be  large  and  expand  u,v,w ,  and  p 
appearing  in  (2.7)  in  the  form 

u  =  uq  4  £  {uj  4  V\C\)  4  0(—\/k 2), 
v  =  kV\C\  4  vq  4  V2C2  4  0(\/k), 

(2.8 a,  6,  c,  d ) 

w  =  W\S]  4  0(1/ A;), 

p  =  k4p0  -f-  fc3p!  +  Pp2  4  k{p\  4  Ci  Pi  }H - • 

Here  Cn  =  cosnA:^,  Sn  =  sin  nkz  and  the  functions  uo,Uo,  etc.  depend  only  on  x  and  y. 
The  Rayleigh  wave  velocity  components  U,  V ,  W  appearing  in  (2.7)  may  be  found  from  the 
large  k  solution  of  (2.5).  In  order  to  determine  this  structure  we  expand  the  wavespeed  c  in 
the  form 

c  =  Co  4  Ci/k  4  (2.9) 

and  a  similar  expansion  is  written  down  for  a.  We  assume  that  u0  ^  <4  in  y\  <  y  <  2/2 •  We 
then  write 

G  =  GkA  4  •  •  •  J=  k\ 

The  appropriate  expansions  for  U,  V,  W  and  P  are 

U  =  tio4}ftC,4-, 

V  =  Vo  4  jVjCi  H - 

(2.10a,  6,  c,  d) 

W  =  pH/i5'i  4  •  •  • 

P  =  Po  4  ■  •  •  4  p-PiCi  4 —  • 

If  we  substitute  the  above  expansions  into  (2.5)  we  find  that  Po  satisfies  the  usual  Rayleigh 
pressure  equation 

Poyy  ~  J^-P0y  ~  <X02 P0  =  0,  (2.1  1«) 

Uq  ~  Cq 


where  Vo  is  given  by 


ia0(u0  -  cq)V0  =  -Poy. 


(2.116) 


At  next  order  we  find  that 


[u0-Co  ( u0-Coy J 


VqU\ 

«0  —  Co 


.j.  2  f*  WOy^lVol 

Wj  =  - -  <  Wly  -  -  \  , 

Wo  -  Co  l.  Wo  -  Co  J 


A  =  2m,  (  VoUly  -  . 

t  Wo  -  Co  J 

We  now  substitute  (2.8)  into  (2.7)  and  equate  the  leading  order  terms  proportional  to  C\  = 
cos  kz  in  the  x  and  y  momentum  equations.  This  yields 


Viw0y  =  -U\ 


\GuaUi  -  =  -V„ 

Wo  —  Co 

These  equations  have  a  consistent  solution  only  if  wo>  Vo  satisfy 

x6g,Bo,-4ltgll,!0|l  =  i- 

Wo  ~  Co 


(2.12) 


This  equation  has  to  be  solved  in  conjunction  with  the  Rayleigh  equation  (2.11a).  In  fact  it 
is  convenient  to  write  (2.11b)  in  terms  of  Pq  using  (2.16)  to  give 


XG(u0)y/2 


4ui 


Oy 


PS, 


<*o(w0  -  co)  (w  -  Co)2 


=  1. 


(2.13) 


In  fact,  since  (2.11a)  does  not  involve  derivatives  with  respect  to  x ,  we  can  multiply  P0  by 
an  arbitrary  function  of  x,  we  therefore  write  Pq  =  ao^f^Po  where  the  eigenfunction  Pq 
satisfies  (2.11a)  together  with 


P oy  =  0,  y  =  0,  P0  -»  0,  y  -»  oo, 


(2.14) 


and  some  normalization  condition.  We  note  also  that  without  any  loss  of  generality  we  can 
take  B,Pq  to  be  real.  We  then  write  (2.13)  as 


xG(ul)y/2 


U0yB2 

Wo  -  Co 


~  2 

p 


Oy 


(w0  -  Co)2 


=  1. 


(2.15) 


Thus  the  outcome  of  the  Rayleigh  wave  forcing  is  to  introduce  the  term  proportional  to  B2 
in  the  local  mean  flow  equation  (2.15).  This  in  effect  means  that  uo  and  P o  can  only  be 
obtained  numerically  because  (2.15)  is  coupled  to  (2.11).  As  in  Hall  and  Lakin  (1988)  we 
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note  that  the  equations  obtained  by  equating  leading  order  vortex  terms  proportional  to  C\ 
do  not  determine  the  vortex.  The  required  equation  is  obtained  by  equating  leading  order 
terms  independent  of  z  in  the  x  momentum  equations,  this  gives 


VqUqx  +  VoXlQ-y  U Oyy  “b  Px  —  o  ^  • 


dy 

This  equation  can  be  integrated  to  give  V,2  in  the  core,  if  the  vortex  is  to  vanish  at  3/],y2  we 
must  have 

/  (u0u0x  -b  VQUQy  -  u0yy  +  Px)dy  —  0. 

Jy  i 

We  shall  now  summarize  the  interaction  equations  and  matching  conditions  which  the  mean 
state  and  Rayleigh  wave  must  satisfy  if  a  vortex-wave  interaction  is  to  take  place. 

For  convenience  we  now  denote  the  zeroth  order  mean  flow  in  Regions  I,  II,  and  III  by 
u(x,y)  =  (u(x,y),v(x,y),  0).  We  also  denote  the  Rayleigh  wavenumber  and  wavespeed  at 
zeroth  order  by  a ,  c  respectively.  The  zeroth  order  problem  for  the  mean  flow  driven  by  the 
vortex-wave  interaction  can  then  be  written  down  in  the  form: 

Regions  I,  II 

ux  “b  Vy  —  0, 

UUX  -b  VUy  =  -p'x  +  Uyy, 


(2.16a,  b ) 


Region  III 


X  “1“  Vy  -  0, 

,  1 

Uy{xGu 

ry2  _ 

/  {itu*  +  VUy  -  Uyy  +  Px }  dy  =  0, 
Jyi 


together  with  the  boundary  conditions 


u  =  v  =  0,  y  =  0, 
u  >  ue,  y  -*  oo, 


and  matching  conditions 


u,  uy ,  u,  P,  Py  continuous  at  y  =  yx ,  2/2- 


(2. 17a, 6) 
(2.18) 

(2.19) 

(2.20) 


We  need  not  write  down  separate  equations  satisfied  by  P  in  I,  II,  III  since,  using  the  notation 
used  above,  the  single  equation 


(u  -  c)(Pyy  -  a2P)  -  2 UyPy  =  0, 


(2.21) 


together  with 

P  =  1,  y  =  2/c,  Py  =  0,  y  =  0,  P  =  0,  y  =  oo,  (2.22) 

where  yc  is  the  critical  layer  location  are  sufficient  to  determine  P.  Finally  we  note  that 
a(i),c(x)  the  eigenvalues  of  (2.21),  (2.22)  must  satisfy 


ac  =  constant, 
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so  that  the  wave  frequency  stays  fixed  as  it  moves  downstream.  In  effect  this  constraint  is 

satisfied  by  adjusting  the  wave  amplitude  B(x)  appropriately  as  x  varies.  We  further  note 

_ 2 

that  thin  shear  layers  of  thickness  k~  are  needed  near  y\,  y2  in  order  to  allow  the  vortices  to 
decay  to  zero  exponentially.  These  layers  are  passive  and  require  only  minor  modifications 
to  the  structure  given  by  Hall  and  Lakin  (1988).  An  analysis  in  these  shear  layers  shows 
that  below  y\  and  above  ?/2  the  vortices  are  smaller  by  a  factor  e~k  raised  to  some  positive 
power.  This  means  that  the  criical  layer  forcing  mechanism  discussed  by  HS  is  exponentially 
small  compared  to  the  distributed  mechanism  discussed  here. 

Before  investigating  the  growth  of  small  amplitude  Rayleigh  waves  from  some  position 
along  the  flow  direction  we  first  discuss  briefly  how  the  eigenrelation  associated  with  (2.21), 
(2.22)  can  be  calculated  when  u  changes  by  a  small  amount.  We  shall  assume  then  that  u 
can  be  written  as 

u  =  u(y)  +  u(y), 

and  we  suppose  u"  is  discontinuous  at  y^.  This  discontinuity  is  associated  with  the  transition 
layers  at  In  these  layers  the  jump  in  u"  is  smoothed  out  in  the  manner  discussed  by 

Hall  and  Lakin  (1988).  However  for  simplicity  we  will  assume  below  that  there  is  just  one 
such  location  in  the  flow,  later  we  will  simply  sum  over  all  such  locations. 

We  perturb  a,  c,  and  P  by  writing 


a  =  a  -f  a, 
c  =  c  +  c, 

yd  =  Vd  +  yd, 

p  =  p  +  p. 


(2.23a,  6,  c) 


We  assume  that  the  quantities  denoted  by  a  tilde  are  smaller  than  those  with  a  bar.  The 
equations  satisfied  by  P,  P  are 


(u  —  c)(P"  —  a2P)  —  2u'P'  =  0, 
Py  =  0,  y  =  0,  P  ->•  0,  y  -*  oo, 


(2.24) 


(tZ  -  c)(P"  -  a2P)  -  2 uP'  ={u-  c){P"  -  a2 P)  +  2 aa(u  -  c)P} 


Py-  o,  y  =  0,  P  -►  0,  y  -»  oo. 


(2.25) 


The  system  (2.24)  constitutes  an  eigenvalue  problem  for  a  with  c  =  u(yc )  with  yc  defined  by 
u"(yc)  =  0-  The  second  system  only  has  a  solution  if  the  appropriate  orthogonality  condition 
is  satisfied.  However  we  confine  our  attention  to  constant  frequency  neutral  modes  so  that 
c  is  given  by 

—u1  u" 

c=--z^  +  uc  (2.26) 


Here  a  subscript  c  denotes  a  quantity  evaluated  at  the  critical  layer.  The  wavenumber 
perturbation  a  is  then  determined  by 


Perhaps  the  easiest  way  to  determine  d  is  to  integrate  the  differential  equation  for  P  once 
by  writing  P  —  PF(y).  In  order  to  write  down  the  condition  that  the  solution  obtained  in 
this  way  is  continuous  across  y  —  yc  we  first  assume  that  in  the  neighborhood  of  yc,  with 


y  =  V  ~Vc 


u  —  <f>o  +  y<f>}  +  y1(f>2  +  •  •  ■ 

-2 

u  =  c+yyi  +  \n-2  +  ■■■■,  /x-2  =  0, 


(2.27a, b) 


p=  ,  ,  y2Q2  fQ* 

2!  3! 


Qi  ==  -a2- 


Here  we  have  normalized  P  to  unity  at  the  critical  layer.  After  some  manipulation  we  find 
that  the  condition  to  determine  d  is 


f<xpS[2aa(n  -c)P  +  2 u'P'  -  ^5^'} 

Jo  ^  (u  —  c)3 


C  DU 

-  TT  -  T 2\dy 

y3  y 1 


-CD 

—j  H - f  E 

iy2c  yc 


(2.28) 


where  C,  D,  E  are  defined  by 


C  =  2Q2fi,  3[c  (j)0] 

D  =  /xr3  [Q3(c  ~  to)  -  • 

E  =  y}PdP^[u"(yZ)  -  u"(yj)}[u{yd)  -  c]' 


(2.29 a,  b) 


Here  we  have  denoted  quantities  evaluated  at  yd  by  a  d  subscript  and  [u"{yd)  —  u"(yj)] 
denotes  the  jump  in  u"  across  yd.  For  a  given  profile  u  we  can  compute  d  from  (2.28)  with 
c  given  by  (2.26),  however  in  ihe  vortex- wave  interaction  we  have  to  maintain  a  constant 
frequency  for  the  wave  as  it  moves  downstream.  Thus  if  u  represents  the  change  in  u  over 
some  small  distance  in  x  then  a  determined  in  this  way  must  satisfy 


dc  +  oc  =  0 


(2.30) 


and  in  general  d,  c  determined  by  (2.26),  (2.28)  will  not  satisfy  this  condition.  Thus  the  role 
of  the  wave  forcing  the  mean  state  now  becomes  apparent  since  B(x ),  the  Rayleigh  wave 
amplitude,  must  now  adjust  in  order  to  enable  us  to  satisfy  this  condition  (2.30).  Now  we 
shall  investigate  the  possible  growth  of  small  but  finite  amplitude  Rayleigh  waves  from  a 
strongly  nonlinear  vortex  flow. 


3.  Small  Amplitude  Rayleigh  Waves  Bifurcating  from  Strongly  Nonlinear  Vortex 
Flows 

In  the  first  instance  we  shall  seek  finite  amplitude  solutions  of  the  interaction  equations 
(2.16)  -  (2.22)  appropriate  to  the  situation  when  the  Rayleigh  wave  is  small.  The  expansions 
we  develop  are  related  to  those  given  by  Hall  and  Lakin  (1988)  for  the  case  when  no  Rayleigh 
forcing  occurs.  The  major  difference  between  the  two  expansion  procedures  is  that  the  initial 
vortex  in  the  work  of  Hall  and  Lakin  was  calculated  from  the  limiting  form  of  the  weakly 
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nonlinear  theory  of  Hall  (1982b).  Before  discussing  how  we  must  incorporate  the  Rayleigh 
wave  into  the  expansions  procedure  we  comment  on  the  solution  of  (2.16)  -  (2.22)  when 
B  =  0. 

Suppose  then  that  at  the  initial  location,  taken  without  loss  of  generality  to  be  x  =  1, 
we  have 

u  =  u{y). 

In  the  absence  of  any  vortex  activity  the  mean  flow  function  u  can  be  expanded  in  powers  of 
A  =  x  —  l,  we  must  however  assume  that  ii//,(0)  =  0  in  order  that  no  logarithmic  terms  occur 
in  the  expansions.  At  first  sight  we  would  expect  that  the  expansions  of  yx  and  y 2  should 
proceed  in  powers  of  A'.  However,  unless  certain  quite  severe  restrictions  are  imposed  on  u, 
it  can  be  shown  that  the  required  continuity  and  boundary  conditions  cannot  be  satisfied. 
This  occurs  because  the  existence  of  a  vortex  A  =  0  means  that  u ,  uy  are  initially  continuous 
at  i/i,i/2.  Thus  the  contributions  to  the  matching  conditions  on  u  a  small  distance  beyond 
X  =  0  arising  from  the  perturbations  to  y\,y-2  vanish  at  0(A).  In  effect  we  must  therefore 
allow  yx  and  1/2  to  be  perturbed  by  0(A!/2)  so  that  the  second  order,  ie  0(A),  correction 
terms  come  into  play.  We  therefore  write 


J/i  —  yi  +  Al/2i/ii  +  Ay  i2  +  •  •  • , 
Vi  —  J/2  +  A 1/2 1/21  +  Xy-ii  +  ■  •  • , 


(3.1a, b) 


and  note  that  the  special  case  yx  —  y2  is  discussed  in  the  next  section.  We  anticipate  a  similar 
expansion  for  the  Rayleigh  disturbance  amplitude  whilst  the  wavenumber  and  wavespeed 
expand  in  whole  powers  of  A".  We  therefore  write 

B  =  X}/2B  +  X Bx  H - , 


a  =  a  +  Ad  +  **-,  (3.2a,  6,  c) 

c  ~  c  +  Ac  +  •  •  • , 

with 

ac  +  ac  =  0,  (3.3) 

if  the  Rayleigh  wave  is  to  remain  neutral  and  of  constant  frequency  as  it  moves  downstream. 
In  Regions  I,  III  we  write 


U  =  U  B 
V  =  VB 


whilst  the  wall  curvature  expands  as 


u  +  X U\  T  A2U2  T  ■  ■  ■ , 
v  +  X  i>\  +  A2i>2  +  •  •  •  , 
7Tq  +  A  7T 1  +  A  2  7T 2  +  •  '  '  . 


\  —  \  0  +  A  \  1  +  A 2  \  2  +  •  •  ■ 


(3.4a,  6,  c) 
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If  we  substitute  the  above  expansions  into  (2.16)  and  equate  terms  of  0(X°)  we  obtain 


U  \+Vy  =  0, 

=  -To  +  Uyy, 


(3.5a,  b) 


so  that  in  Regions  I,  III  the  required  solutions  are 

u i  =  -  u'  j  +  J  (7r0  -  u")u-2dy j  , 


v  —  JjU  +  u  [tto  —  u"\u  2dy, 
Jy, 


with  j  =  1,2  and  Jj,  J2  constants  to  be  determined. 

In  fact  Ji  is  fixed  by  the  condition  u\  =  0,  y  =  0,  this  gives 


(3.6a,  b) 


fVl 

Ji=  (tto  -  u")u~2dy.  (3.7) 

Jo 

and  we  note  that  the  above  integral  exists  if  we  assume  that  u"  —  7r0,  u'"  =  0,  y  =  0.  In  the 
core  region  III  we  retain  (3.4)  and  expand  the  Rayleigh  pressure  function  P  as 


P  =  P  +  XPt  +  ■  ■  ■  . 


At  zeroth  order  the  coreflow  downstream  velocity  component  is  given  by 

\/A  -f  2  y 

and  the  zeroth  order  approximations  to  (2.17a,b)  yield 


1  (  rs  ni  n  Xiv/ A  +  2 y 

1  _  1A  B  2xo\/xoCj  ’ 


d  =  b-  ^  - -■ N- -  + 

/w-(xoGp/; 


B2  p 


where 


Xi(A  +  2y)2 

6Xo(XoG')^ 


TW  -  c(xoG)V2)  •  (3.8a,  6,  c) 

Here  the  constants  A,  B ,  K  and  b  are  to  be  determined.  In  fact  the  zeroth  order  approxima¬ 
tions  to  the  condition  that  v  is  continuous  at  yj  yield 


y/A+  2t/i(Ji  +  K) 


-  by/Uxo  +  6^(A  +  2j/i)5  , 


v/A  +  2j/2(  J2  +  K  —  B2 12) 


bJo^o  -  B2  r  jA  +  2yI'dy  +  -^-(A  +  2y,)l. 

''yi  OXo 

(3.9 a,  6) 
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In  addition  the  continuity  of  u,uy  at  zeroth  order  at  y  =  yi,j/2  yields 

\/ A  2 ~y~j  1 


(3.10a, 6) 


We  recall  that  </]  must  satisfy  (3.7)  in  order  that  wi(0)  =  0,  then  (3.9a,b)  determine  b  and 
J2  in  terms  of  one  remaining  unknown  constant  I\ .  We  now  use  the  coreflow  expansion  to 
show  that  the  zeroth  order  approximation  to  the  vortex  condition  (2.18)  is 


r~c  b^A  +  2y  +  1  iA  +  2y)2x  1  1 

VxoG  12  XqG  \/ A  +  2y%/xoG 


5-Jro  y  =  B2I{S2-S ,)■ 


where 


S]  =  JA  +  2yj. 


I(S2  -  5,)  =  jT  I(y)(\/A  +  2y2/y/A  +  2y  -  2)dy. 


Hence  the  above  equation  together  with  (3.7),  (3.9a,b)  yield 


b  =  6,  -  £2/, 


4'  - 


ts  _  WxoG 

K  =  ~sT~ 


kqG  ry'  — 2.  ,  Xi^i 

7 - /  (tt0  -  w  )u  rfy  +  -7—, 

>1  7o  OXo 


(3.11a,b,c) 


,,  ,  & Sg  VX+lyl'dy  XI 53 

J2  =  —A  +  £>  /2  H - 7, - 7, - T  ~  ■ 

62  C>2  OXo 


We  now  substitute  for  itj  from  (3.6a)  into  (3.4a)  and  expand  about  y  =  t/1,1/2  to  obtain 
u  =  u3  +  X1/2u0yjt  +  X  |u°j/j2  +  u"3^—  +  — -ry — -  -  uj Jj|  4 - , 

*•  =  u'3  4-  A’I/2u"J'y,-1  +  X  { u"'fc2  4-  u"'%  A  ( U"-~  — ^  -  - 


Uy  =  u3  4-  Xl'2u"3yji  A  x  |u'0j/j2  4-  u'"3 4-  ^ — — — J  -  — -  «'°]  -  Jju 4-  •  •  •  . 

(3.12a,  b) 

Here  u'\  u"\  etc.  denote  u'(ijj),  u"(2/j),  etc.  and  the  next  order  correction  terms  are  0(X3^2). 
If  we  perform  similar  expansions  for  the  coreflow  solution  we  obtain 


u  - 


Uy  = 


S,  ,  X'^j,  ,  y  f  1  ....  ,  -  6,  .  |  XxSj _ ! £ill, 

x&  +  Vx&s<  +  *\^asr+m  B,jl  2x„v^c  :\  x, 


V  XoG  VxoG  Sj 

1  X'!%y 

V  Xo  GSj  \JxoGS3} 


3  ■* - , 


(3.12c,  d) 
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If  the  initial  profile  is  chosen  such  that  u ?  =  Sj/ y/xoG,  u'J  —  1  jSj  \J\qG)  we  see  that  the 
order  X°,  X1 /2,  terms  in  (3.12a)  and  (3.12c)  are  consistent  whilst  also  are  the  order  X0  terms 
in  (3.12b),  (3. 1 2d).  However  if  we  are  to  make  u  and  uy  consistent  up  to  orders  X,  X1/2 
respectively  then  we  have  four  equations  to  satisfy  but  only  three  unknowns  yu,  y2\,  B2  at 
our  disposal.  Thus  we  cannot  in  general  make  u,uy  consistent  at  these  orders  and  the  only 
possible  remedy  is  to  allow  for  diffusion  layers  of  depth  Xx^2  in  the  neighborhood  of  y\,y2. 
However  these  layers  occur  only  in  I,  II  where  the  mean  flow  satisfies  the  boundary  layer 
equations.  These  blending  layers  occur  where  d2  ~  dx  and  enable  jumps  in  u,uy  across  the 
layers  to  be  accommodated. 

We  shall  only  consider  the  behavior  near  y  =  y2,  a  similar  structure  holds  at  y  =  y\.  We 
define  a  variable  £  by 

so  that  We  modify  (3.4a,b)  near  y  =  y2  by  writing 


so  that  ~ 

OX 


u  =  uB  +  XU(0- +  •••, 

v  =  +  X1'2^)  +  •  •  • , 

and  it  is  straightforward  to  show  that  U ,  V  satisfy 

U"  =  2  ((/-£{/'), 

(rfb) 72  V'  =  (2  ~  #)  • 

The  solution  of  (3.14c)  which  decays  when  £  — >  00  is 

tie-*14  =  C2e-?I4U  (|{)  , 


(3.13a,  b) 


(3.14c,  d) 


(3.15) 


where  U (a,  x )  denotes  the  parabolic  cylinder  function,  C2  is  a  constant  and  £>  =  y2 1  ( -Qxab ) 
The  effect  of  this  extra  term  in  the  expansion  of  u,  uy  near  y2  is  to  produce  terms 

C2XU  (§>£2)  ,  —  C2Xx I2  [2Jxa&)  ^  ^  (l’^2)  on  r*8ht  hand  sides  of  (3.12a,b)  for  j  =  2. 
A  similar  analyses  near  y  =  y  1  produces  similar  terms.  We  are  now  in  a  position  to  make 
(3.11a),  (3.12a)  consistent  up  to  order  X  and  (3.11b),  (3.12b)  consistent  up  to  order  X1^2. 
This  is  achieved  if 


{“'"f  +  +  c’u (§’(~,)%)  “ 


2xo\/xo^ 


— 3/ji 

y/xo^Sj' 


(3.16a, b) 
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If  CUC2  are  eliminated  above  we  obtain 


=  OTi1*  -  ^  -  iikz  ~  +  ^  <3l7<l) 


'  2  +2^usr  ^-2{  l)u(u-m,)  fi-,  ■ 


(3.176,  c) 


In  the  absence  of  a  Rayleigh  wave,  ie  when  B  =  0,  we  can  substitute  for  K  from  (3.12b)  with 
b  given  by  (3.12a)  to  give  a  quadratic  equation  for  y:\.  In  the  presence  of  a  Rayleigh  wave 
we  can  only  express  yj\  in  terms  of  B 2 .  However  if  we  now  derive  the  condition  that  the 
wave  remains  neutral  up  to  O(X)  then  a  further  equation  linking  y:\,B2  will  be  obtained. 
In  addition  we  note  that  if  the  initial  data  is  such  that  Aj  vanishes  for  j  =  1,2  then  the 
expansions  for  yj  proceed  in  whole  powers  of  X. 

If  we  use  the  result  (2.29)  applied  to  the  O(X)  correction  to  the  mean  velocity  field  we 
find  that 


fc,  +  Khi  +  B2h3  +  J,h,  +C,h„U(2;,  (-1)'6)  +  C2heU(-,{-\y(2)  =  0. 


(3.18) 


Here  the  constants  h\y  h2,  h3 ,  h4,  hs,  h6  are  defined  by 


/•<*>  1  =  f  2 a&P  2 u'cP'  \ 

Jo  [  \  (u  —  c)2  (it  —  c)4  J 

Jo  {u  —  c)3  [  u  —  c  J 


(y  -  Vcf  ( y  -  yc)2 


2  y2c  Vc 


ry-i  PP'  (c  -  2 u)  , 

dy  +  X  i/  7r - rrj- — jTTT-dy 

Jii 1  (u  -  c)4  Xo^u 


+yiAP?)'ki [u(yi)  -  c]  3  -  y2i=(f,22)'A2My2)  -  c]‘ 
r  2PP'  (  u'M  \  , 

+L  Ffr 


m  2rr  f  ^.c  rtl  . 

WSl  (rrE)?b»G“} 

rv «  2PF  f  7[c  —  2n]  /'  1 

3  Jy ,  (u  -  c)3  \  u3XoG2(u  -c)  +  xoGu  J 

,4  =  _r 

7y2  [U  -  C)3  (  U  -  C  J 


2PiP{ 

(u(y,)  -  c)3’ 

2P2P± 

(«(&)  -  c)3’ 
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with 


M(y)  = 

M(y)  = 


u"  -  7T0  fy  (  7Tq  -  u"\ 

u  Jo  \  u2  J 

u"  -  7T0  fy  ( 7T0  u"\ 

— — uL  {-*-) 


dy , 


dy. 


and  C,  D  given  by  (2.29)  with  c—u\c—  and  u,  given  by  (3.6a).  Furthermore  Pj,P j 
etc.  denote  P,P'  evaluated  at  yj.  If  we  eliminate  C\,C2  and  J 2  from  (3.18)  using  (3.16b) 
and  (3.12c)  we  obtain 

/17B2  =  hg  +  h5u>\A.iyu  -f  h^L02^2V2\  (3.19) 


with 


/l7  — 

hs  = 


»  yJA  +  2 yl’dy 


+ |J + /.3 + /.« |/2  -  /ti 
-K  -  k  _  7,}  -  + 


X\SL 
6x0  J 


/l2  —  /l2  —  /I4. 


Finally  we  then  substitute  for  B 2  from  (3.19)  into  (3.17a)  with  j  =  1,2  to  give  the  coupled 
equations  for  yu,yi\‘ 


A2  jj/jl  +  £21^2  ~  | 


hio  +  hghah7^  +  uj2A.2h6h.gh1 1  j/21 » 


A12  4"  hiihghj^  +  Ajhshjjhj  1y1j , 


(3.2 Oa,  b ) 


with 


^>1 


&i_  _  xi _  («"*  ~  ftp) 

10  52  ZxoVXoG  u1 


/  /+  ,  61  Xi&  (u"2-7r0) 

S22  v/x^ST  12  S22  3xon/x^  «» 


(3.21a,  6,  c,  d) 


and  7+defined  by 

/+  =  /V2  I'yjA  +  2ydy 
Jy  1 

Thus  the  initial  small  amplitude  form  of  the  vortex  and  wave  are  determined  at  leading 
order  when  the  nonlinear  equations  (3.20a,b)  are  solved  for  yn,y2 1-  Note  however  that  /j8 
appearing  above  is  a  quadratic  function  of  these  quantities  whilst  u>uu>2  depend  on  ynj/21, 
through  parabolic  cylinder  functions  so  that  the  solution  of  (3.20)  must  be  found  numerically. 
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A  numerical  example 

In  order  to  proceed  further  we  must  specify  an  initial  velocity  profile  and  curvature  distribu¬ 
tion.  Without  any  loss  of  generality  we  may  take  Xo  =  1  and  we  then  take  the  initial  profile 
u  to  be  given  by 


«  =  +  ^(y2  -  ^-)J  ,  0  <  y  <  y,  =  1 


u  = 


\J  A +  2y 


Vg 


y\  <  y  <  y*  =  2, 


\JA  +  ^h  \ 

u=-7=-+1-exp| 


y-y 


\ 


y2  <  y  <  oo 


(3.22a, i>, c) 


\J  A  +  2y2\fG  j 

Here  the  constant  A  is  to  be  specified  but  we  restrict  our  attention  to  situations  when  u'  is 
always  positive.  In  order  that  u  and  uy  are  continuous  at  the  positions  y  =  yi,y2  we  require 
that  A  and  G  are  defined  by 


A  = 

G  = 


~~(A  +  e) 

(A  +  i)  ’ 

1 


(A  +  IKA  +  I) 


The  profile  defined  above  has  an  inflection  point  at  y  =  \  so  that  the  wavespeed  is  given  by 


__  A 

C~  2  +  192' 


(3.23) 


The  first  step  in  the  calculations  is  to  solve  Rayleigh’s  equation  for  P  with  u  and  c  as  given 
by  (3.22)  and  (3.23).  In  Figure  (3.1)  we  show  the  wavenumber  a  as  a  function  of  A.  The 
eigenfunctions  P  associated  with  these  modes  are  shown  in  Figure  (3.2).  Though  it  cannot 
be  seen  easily  in  this  figure  we  note  that  at  the  transition  layers  yi,y2  the  quantities  P,  P' 
are  continuous  whilst  P"  is  discontinuous. 

The  constants  hi,h2,  etc  appearing  in  (3.21)  can  then  be  determined  using  Simpsons  Rule 
to  evaluate  the  integral  and  a  routine  to  evaluate  parabolic  cylinder  functions.  This  was  done 
using  either  the  series  or  large  argument  asymptotic  form  of  the  parabolic  cylinder  functions 
dependent  on  the  size  of  the  argument.  The  equations  (3.21)  were  then  solved  using  a  Newton 
iteration  procedure.  The  calculations  were  restricted  to  the  case  A  =  .5,1, 1.5, 2, 2.5,3.  but 
different  values  of  the  curvature  parameter  \i  were  used.  In  Figures  (3.3a,b,c)  we  show 
the  computed  values  of  jq,  y2  and  B 2  as  functions  of  xi-  The  calculations  were  carried 
out  for  the  values  of  A  used  to  calculate  Figures  (3. 1,3. 2)  but  in  fact  we  were  unable  to 
find  solutions  of  (3.21)  for  the  case  A  =  .5.  Furthermore  the  solutions  shown  could  not  be 
extended  beyond  the  ranges  of  Xi  shown.  At  the  lower  end  of  the  ranges  in  question  the 
results  showed  that  yu,y 2i  were  tending  to  a  constant  plus  a  multiple  of  the  square  root  of 
the  difference  between  Xi  an<i  it’s  cutoff  value.  This  suggested  to  us  that  at  such  a  point 
two  solutions  of  (3.21)  were  coalescing,  however  careful  searches  were  unable  to  find  a  second 
root.  However  it  could  well  be  that  other  roots  exist  but  our  Newton  iteration  procedure 
was  not  able  to  detect  them.  Since  the  main  aim  of  the  present  calculation  was  to  determine 
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whether  physically  acceptable  solutions  of  the  vortex-wave  interaction  equations  could  be 
found  near  the  position  of  neutral  instability  of  a  Rayleigh  wave  we  did  not  investigate  the 
matter  further. 

In  fact  an  inspection  of  Figure  (3.3c)  shows  that  B 2  is  negative  for  a  range  of  values  of 
These  solutions  are  not  physically  acceptable  since  B  must  be  real.  At  first  sight  it  might 
appear  that  these  solutions  are  relevant  to  the  case  when  the  Rayleigh  wave  is  bifurcating 
subcritically  from  x  =  1.  This  is  not  the  case  since  our  expansion  procedure  cannot  be 
simply  modified  to  take  account  of  x  —  1  being  negative.  This  is  because  the  blending  layer 
structure  now  has  solutions  increasing  exponentially  away  from  y\,yi  so  that  the  matching 
of  u,uy,  v  cannot  be  achieved  where  /, //  and  //,///  meet.  Nevertheless  we  feel  that  the 
solutions  of  (3,21)  with  negative  B 2  do  have  some  physical  significance;  we  return  to  that 
point  in  the  final  section  of  this  paper. 


4.  The  Spontaneous  Generation  of  Rayleigh  Waves/ Vortices 


Suppose  that  we  have  an  incoming  inflectional  velocity  profile  u  =  ii(y)  at  .r  =  1.  We 
use  the  notation  of  the  previous  section  and  denote  the  critical  level  by  y  —  yc  and  denote 
quantities  evaluated  at  the  critical  layer  by  a  subscript  c.  We  suppose  further  that  the 

position  y  =  y  at  which  Guuy  —  1,  where  finite  amplitude  vortices  emerge,  is,  without  any 

loss  of  generality,  such  that  y  >  Vc  >  1. 

In  the  first  instance  we  ignore  the  mean  flow  correction  driven  by  the  vortex  structure 
and  write 

tx  —  u  X  -F  X2U2  T  •  •  • , 

v  =  v  +  Xu,  F  X^v-2  +  •  •  • ,  (4. In,  6,  c) 

P  =  7To  +  X7Tj  T  X27T2  +  •  •  •  , 

where  X  =  x  —  1.  The  expansions  for  tx,  v  fail  when  y  =  0{Xy^)  where  a  passive  boundary 
layer  is  needed  to  satisfy  the  no-slip  condition.  It  should  be  noted  here  that  v  cannot  be 
specified  arbitrarily  but  must  be  determined  in  terms  of  the  initial  streamwise  profile  u{y). 
If  the  above  expansions  are  substituted  into  the  boundary  layer  equations  then  we  can  show 
that 

«i  +  v'  =  0,  tux,  +  vu  =  —  7T0  +  u",  (4.2a,  6) 

and 

2tx2  +  v\  =0,  2tux2  -F  txj  +  vu\  +  v\u  =  —7 r,  -f  u".  (4.3n,  b ) 

These  equations  can  be  integrated  to  give 


ry 

v  =  tx  /  [7T0  -  u"]u~2dy, 

Jo 

ry 

tx,  =  —11  I  [7T0  —  tx"]tx-2<x'y  + 
Jo 


(4.4a, b) 


tto  —  u 


u 
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and 


fV 

t}j  =  U  /  [ 7T!  —  u "  +  u2  +  vu\]u~2dy, 

Jo 

.  fyr  -2  -//  _ /  -i _ 2  i  ,  {7Tl  -  u'[  +  U2  +  Vu\} 

U2  =  —U  [tt,  +  u\  -  ux  +  vu\)u  2dy  +  - - r - 

Jo  u 


(4.4a,  b) 


Here  we  have  made  assumptions  about  the  behavior  of  u  near  j/  =  0  in  order  that  the  integrals 
shown  exist  at  y  =  0.  If  these  assumptions  are  not  made  then  the  solution  for  y  =  0(1)  must 
have  terms  involving  log  X  in  its  small  X  form.  In  view  of  the  discussion  given  at  the  end 
of  the  section  2  it  follows  that  a  and  c  for  the  Rayleigh  wave  associated  with  the  perturbed 
velocity  profile  (4.1)  must  expand  as 


a  =  a  +  di  X  +  a2X2  +  -  •  • , 
c  =  c  +  c-iX  +  c2X2  4 - . 


(4.5a,  b) 


If  we  use  the  notation  of  §2  then,  in  the  absence  of  vortex  activity,  it  follows  that  &i,  c\ ,  <*2,  c2 
are  determined  by 


r 

Jo 


-  a)P  +  2fi'i^P-  +  2u'P'  -  cf  -  d  =_CL+Di 

y  2  y2c  yc 


(u  —  c)2 


J  =  1,2.  (4.6) 


Here  C:,  D3  satisfy  (2.29)  with  (c  -  4>o)  replaced  by  [cj  —  uj(yc)).  However  if  the  wave  is  to 
remain  neutral  and  of  constant  frequency  we  also  require  that 


Cl 


c2 


-u'cu" 


1 C 


u'c 


+  Ulc 


*'c 

\  <  2c 


u 


III 


+  u2ci 


(4.7a,  b,  c,  d) 


a^  +  acx  =  0, 


aiCi  +  a2c  +  ctc2  =  0. 


The  conditions  (4.6),  (4.7)  impose  a  constraint  on  the  initial  profile  u  =  u(y)  if  the  Rayleigh 
wave  is  to  stay  neutral  up  to  order  X2.  Shortly  we  shall  see  how  a  small  amplitude  vortex 
emanating  from  y  =  y  changes  this  conclusion.  An  examination  of  the  passive  boundary 
layer  where  y  =  0(X 1|/3)  shows  that  this  layer  does  not  influence  (4.6),  (4.7). 

Now  let  us  find  the  effect  of  a  finite  amplitude  vortex  on  the  above  discussion.  We 
recall  from  the  work  of  Hall  and  Lakin  (1988)  that,  if  Guuy  =  1  at  x  =  1,  then  a  finite 
amplitude  vortex  will  initially  grow  within  the  region  —  XA^2y\  +  y  <  y  <  Xl^2yi  +  y  for 
small  values  of  X.  A  result  of  some  consequence  for  the  present  investigation  is  that  the 
mean  flow  modification  driven  by  the  vortex  is  0(X3^2)  and  confined  to  the  same  region 
as  the  vortex.  We  shall  choose  B,  the  Rayleigh  wave  vortex  amplitude  to  be  of  that  size 
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which  also  leads  to  a  mean  flow  modification  of  size  X 3/*2.  The  corresponding  effect  of  this 
mean  flow  correction  on  a  and  c  can  then  be  again  deduced  from  (2.29)  and  (2.30).  In  fact 
the  size  of  this  effect  can  be  found  directly  from  (2.29)  by  noting  that  u  is  now  a  function 
confined  to  a  layer  of  depth  0{XX^2)  well  away  from  the  critical  layer.  This  means  that  the 
only  contribution  to  (2.29)  from  u  comes  from  the  terms  in  the  integrand  proportional  to 
P1 .  Since  u  ~  0(X 3//2)  we  should  anticipate  a  contribution  of  this  size  to  (2.29),  however 
we  shall  see  below  that  some  cancellation  in  the  integral  occurs  so  that  only  a  contribution 
of  size  0(X s/2)  is  generated.  Moreover  since  u  is  exponentially  small  away  from  the  critical 
layer,  the  order  X5^2  correction  terms  in  the  expansions  of  a  and  c  must  be  identically  zero. 
In  this  case  the  satisfaction  of  (2.29)  is  achieved  by  the  adjustment  of  B(X). 

We  define  a  variable  9  by 

(4-13) 


where 


Gu  uy  -  1,  (u  uy)y  =  0,  y  =  y. 


We  shall  seek  a  solution  of  (2.16)  -  (2.20)  by  expanding 

B  =  BXxI2  +  --- 


(4.14) 


and  the  velocity  field  in  the  form 

u  —  Ug o  +  X1/20uio  +  X  {02U2o  +  u0i}  +  Xz!29u\\  +  X3^2u\](0)  4-  •  •  ■ , 

V  =  Vqq  -f  Xl/20vio  +  X  {92v  20  -f  Uoi}  +  X3/20v)  i  +  Xvm(0 )  +  •  •  • ,  . 
Here  the  constants  uoo,  etc.  are  defined  by 

uoo  =  u(y),  uio  =  u'(y )  =  a~—,  u20  =  2 u"(y), 


(4.15a,  b) 


uoi  =  ui(y),  uu  =  u\(y), 


uoo  =  v{y),  u10  =  v'(y), 


(4.16a  —  i) 


V20  =  2v"(y),  v0i(y)  =  ^(y),  vu  =  uj(y). 

Furthermore  we  assume  that  the  wall  curvature  X  expands  locally  as 

X  =  Xo  +  Xx\  +  X2X2  +  ••••  (4.17) 

The  Region  III  in  this  new  notation  is  defined  by 

-0X1/2  <9  <  9Xx/ 2 

whilst  I,  III  are  determined  by  9  >  9Xx!2,  9  <  —9XxB  respectively.  Here  9  is  a  constant  to 
be  determined  and,  because  of  symmetry,  it  will  be  sufficient  for  us  to  consider  only  9  >  0. 
In  Region  II  we  find,  by  equating  successive  powers  of  X1/2  in  (2.17a),  that 

Gx  ouoouio  =  1, 

UjQ  +  2«0OW2O  =  0, 
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—  Xl  9um  302Uio«2O  .  1*01^10  f)2 

77-^ —  =  -57-  +  Wll  + - + - 07^ 

(jXqUqo  ov  Woo  Woo 


where 


7  =  7  “10  -1  {P'\u  ~  c)  2}  ,  (4.18«,6,  c,  d) 

[Woo  —  CJ  L  J 

evaluated  at  y  =  y.  Equations  (4.18a,b)  are  automatically  satisfied  because  of  (4.14)  and 
(4.18c)  can  be  integrated  once  to  give 


UM  —  {Ao  +  1  B2}9  —  A1#3, 


(4.19) 


X'lUio  U01W10  .  W10W20 

Aq  —  — wu - ,  Aj  =  - . 


(4.20) 


We  have  assumed  above  that  um  is  an  odd  function  of  9  as  was  found  to  be  the  case  by  Hall 
and  Lakin  (1988).  Following  Hall  and  Lakin  we  find  that  in  Region  I 


um  = 


{X0e  +  7 B26  -  A,0‘3)exp  {  — ^[02  -  02]}  U  (1,  yffio) 
=  <'(!.,/¥*) 


(4.21) 


where  we  have  already  made  u  continuous  at  the  junction  of  I, II  and  the  continuity  of  uy  at 
the  junction  of  I,  II  then  yields 


=  -v- 

V  2  \,/(l,  STS)}  »[A0 


uibjr°) 


j  +  7  B2  —  A 162 


(4.22) 


In  the  absence  of  the  Rayleigh  wave  forcing  we  can  set  B  =  0,  equation  (4.22)  then  reduces 
to  the  nonlinear  eigenvalue  problem  for  0  found  by  Hall  and  Lakin  (1988).  Here  equation 
(4.22)  is  not  sufficient  to  determine  B 2  and  0,  the  required  extra  condition  is  obtained  by 
insisting  that  the  0(X3B)  mean  flow  correction  in  I,  II,  III  does  not  make  the  Rayleigh 
wave  wavenumber  or  wavespeed  vary.  The  effect  of  the  mean  flow  correction  can  be  seen 
from  equation  (2.29).  We  expect  that  the  term  in  (2.29)  proportional  to  u'  will  provide  the 
dominant  contribution  to  the  integral,  this  would  suggests  that  the  dominant  contribution 
to  the  integral  would  be  of  order  Xi.  However  this  contribution  vanishes  because  um  itself 
vanishes  at  infinity.  Similarly  the  order  X2  contribution  vanishes  because  um  is  an  odd 
function  of  9.  Now  we  must  set  6t  =  c=A  =  B  =  0  because  the  mean  flow  correction  is 
exponentially  small  at  the  critical  layer.  Thus  (2.29)  leads  to  the  equation 


fOO 

/  0uM{9)d0  =  0. 

J  OO 


(4.23) 


We  note  also  that  the  higher  order  mean  flow  correction  terms  also  vanish  at  infinity  and 
are  odd  in  9  and  therefore  do  not  contribute  to  (4.23).  I11  the  derivation  of  (4.23)  we  have 
used  the  result 

/OO  too 

92u'M(0)d9  =  —2  /  9uMd9 

•OO  J  —OO 
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Equation  (4.23)  yields 


6  f  ,  ,  4  / 11  V^)  1  _  7^2  +  A0  -  |Ai<92 

502uOo  \  V  2  u  (l  y/^B)  J  7^  +  A0  -  \\02  ' 


(4.24) 


Thus  #  and  B  are  determined  by  the  coupled  system  (4.22),  (4.24),  the  system  is  simplified 


by  writing 


(4.25a,  b) 


in  which  case  we  obtain  the  coupled  system 


7b2  +  AQ-3A!©2  =  e^(f,e) 

7£2  +  ao-a102  £/(f,©)  ’ 


7B2  +  Ap  -  |Ai02  _ _ 3_  |  ,0  u  (2^)  ) 

752  +  ao-a,02  502  }  i/(|,e)  /’ 

and  if  we  multiply  (4.26b)  by  5  and  add  to  (4.26a)  we  obtain 


(4.26a,  b) 


—3(1  +  Qa)  _  v{ f.e) 

0(3  +  02)  v(\ ,e)' 


(4.27) 


The  right  hand  side  of  the  above  equation  is  always  positive  so  there  are  no  roots  of  the 
equation  for  0  positive.  Since  we  have  implicitly  assumed  above  that  6  is  positive  we 
conclude  that  there  are  no  acceptable  finite  amplitude  solutions  of  the  interaction  problem. 
Thus  we  have  found  that  small  wavelength  vortices  and  Rayleigh  waves  cannot  spontaneously 
be  generated  in  a  centrifugally  and  inviscidlv  unstable  boundary  layer.  We  conclude  that 
the  only  possibility  is  that  described  in  the  previous  section  where  we  showed  that  small 
amplitude  Rayleigh  waves  can  be  generated  from  an  0(1)  vortex  field. 


5.  Further  discussion  and  conclusion 

In  Section  3  we  saw  that  the  initial  stages  of  vortex-wave  interactions  in  highly  curved 
boundary  layers  can  be  expressed  in  terms  of  an  asymptotic  expansion  in  powers  of  the 
square  root  of  the  distance  from  the  point  of  neutral  stability.  Our  numerical  investigations 
of  a  class  of  initial  profiles  showed  that  in  some  cases  the  interaction  cannot  occur  because 
the  predicted  amplitude  of  the  Rayleigh  wave  is  imaginary.  In  nonlinear  hydrodynamic 
instability  theory  this  would  usually  suggest  that  a  subcritical  finite  amplitude  instability 
occurs.  Here  this  is  not  the  case  because  the  blending  layer  structure  fails  when  x  <  1 
because  the  parabolic  cylinder  functions  now  grow  exponentially  away  from  the  transition 
layers. 

Thus  the  solutions  obtained  in  Section  3  with  B2  negative  cannot  be  used  to  begin  a 
marching  procedure  to  solve  the  full  interaction  equations.  Nevertheless  we  believe  that 
these  solutions  are  still  of  some  interest  and  are  indeed  of  particular  physical  importance. 
This  is  because  these  solutions  will  play  a  crucial  role  when  the  streamwise  lengthscale  x  —  1 


becomes  sufficiently  small  for  nonparallel  effects  to  come  into  play;  see  Hall  and  Smith  (1984) 
for  a  related  discussion  in  connection  with  Tollmien-Schlichting  waves  in  growing  boundary 
layers.  We  do  not  address  the  nonparallel  problem  here  but  we  note  that  nonparallel  effects 
come  into  play  when  x  —  1  becomes  so  small  that  streamwise  derivatives  of  the  Rayleigh 
wave  amplitude  balance  with  changes  in  the  amplitude  induced  by  the  variation  of  the  mean 
state.  In  this  case  the  interaction  is  governed  by  an  integro-differential  equation  whose  large 
x  asymptotic  form  reduces  to  the  small  x  —  1  form  found  in  Section  3.  We  believe  that 
the  negative  B 2  solutions  found  in  Section  3  lead  to  solutions  of  the  nonparallel  evolution 
equations  having  a  finite  distance  singularity.  This  would  explain  why  they  do  not  connect 
with  small  x  —  1  solutions,  whilst  we  believe  that  the  solutions  obtained  in  Section  3  with 
positive  B2  can  be  connected  with  solutions  of  the  evolution  equations. 

In  Section  4  we  investigated  the  special  case  which  occurs  when  the  interaction  begins 
with  the  vortex  and  Rayleigh  wave  both  having  small  amplitudes.  Here  we  found  that  the 
nonlinear  eigenrelation  (4.27)  does  not  have  any  real  solution  so  that  the  interaction  cannot 
take  place.  We  conclude  that  the  spontaneous  generation  of  Rayleigh  waves  and  vortices 
cannot  occur  so  that  the  only  possibility  is  that  discussed  by  Hall  and  Lakin  (1988)  with 
just  a  small  amplitude  vortex  emerging  from  the  neutral  position.  However  there  is  again 
a  nonparallel  evolution  problem  associated  with  the  analysis  of  Section  4  if  X  becomes 
sufficiently  small.  We  expect  that  on  the  appropriately  reduced  streamwise  lengthscale  the 
structure  outlined  in  Section  4  can  be  modified  to  show  that  the  spontaneous  generation  of 
wave  and  vortices  leads  to  a  finite  distance  singularity. 
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